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', Abstract 

>: 

We study the d u asymmetry of the proton in a model approach recently 
Q_i' developed, in which hadronic fluctuations of the nucleon are generated through 

■ gluon splitting and recombination mechanisms. Within this framework, it is 



shown that both d/u and d — u distributions in the proton can be consistently 
described by including only nucleon fluctuations to \ttN) and \irA) bound states. 
Predictions of the model closely agree with the recent experimental data of the 
E866/NuSea Collaboration. 

In 1991, the New Muon Collaboration (NMC) |I] presented a determination of the 
non-singlet structure function i*f — F% at Q 2 = 4 GeV 2 over the range 0.004 < x < 
0.08. From this measurement, the Gottfried Sum Rule (GSR) p| was estimated and a 
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significantly lower value than the 1/3 predicted by the quark-parton model was found. 
Although it could be due to an abnormal behavior of the valence quark distributions in 
the unmeasured region ||, the above result was attributed to a d— u asymmetry in the 
light nucleon sea (see e.g. Ref. 0]). Later on, the NA51 Collaboration [0] determined 
the value d/u = 1.961 ± 0.252 at x — 0.18 in Drell-Yan dimuon production, giving 
strong experimental support to the d — u asymmetry explanation. 

Most recently, the E866/NuSea Collaboration measured the ratio d/u in the nu- 
cleon over the range 0.02 < x < 0.345 in Drell-Yan dimuon production from p — p 
and p — D interactions |J. From this measurement, the E866/NuSea Collabora- 
tion extracted the x dependence of d — u and estimated a value for its integral, 
Jq 1 dx[d — u] = 0.100 ± 0.018, indicating a strong GSR violation. Notably, this value 
is only about two thirds of the NMC estimate for the same integral, 0.147 ± 0.026. 
We will address this issue later on. 

In conclusion, although a d u asymmetry in the nucleon sea has been firmly 
established from experiments, its origin and precise features are yet unclear. Several 
ideas have been put forward to try to explain the GSR violation and the d — u 
asymmetry in nucleons. Among them the Pauli exclusion principle, which would 
inhibit the development of up (down) quarks and anti-quarks in the proton (neutron) 
sea, a pioneer idea by Field and Feynman 0; fluctuations of valence quarks into 
quarks plus massless pions ||, an effect which is calculable in Chiral Field Theory; 
and earlier versions of the Pion Cloud Model (PCM) ||. However, none of these 
attempts gave a satisfactory description of the experimental status of both d — u and 
d/u distributions. It should be noted that previous versions of the PCM used rather 
hard pion distributions inside the nucleon resulting in large u and d distributions 
beyond x ~ 0.3. These large pion contributions can not be easily compensated to 
conveniently describe the fast fall-off of d — u in the whole x range. Albeit large for 
small x, the d — u distribution seems to be negligible beyond x ~ 0.3. In addition, the 
d/u ratio predicted by these models exhibits a dramatic growing behavior not seen 
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in the experimental data. For reviews, see Ref . [|10| 



In this work we will show that a recently proposed version of the Pion Cloud 
Model (PCM) |l]J allows a remarkable prediction of the nucleon's d u asymmetry in 
accordance with the recent results of the E866/NuSea Collaboration. Our approach 
is based on both perturbative and effective degrees of freedom, and it relies on a 
recombination model description of the hadronic fluctuations of the nucleon. 



Let us briefly recall the model introduced in Ref. We start by considering a 
simple picture of the ground state of the proton in the infinite momentum frame as 
formed by three valence quark clusters or valons [|T2|]. The valon distributions in the 
proton are given by 

v(x) = T^v^(l , (1) 

where, for simplicity, we do not distinguish between u and d valons. 

The higher order contributions to the proton structure are identified with meson- 
baryon bound states in an expansion of the nucleon wave- function in terms of hadronic 
Fock states. Such hadronic fluctuations are built up by allowing that a valon emits a 
gluon which, before interacting with the remaining valons, decays perturbatively into 
a qq pair. This quark anti-quark pair subsequently recombines with the valons so as 
to form a meson-baryon bound state. 

The probability distributions of the initial perturbative qq pair can be calculated 
by means of the Altarelli-Parisi Jl3fl splitting functions: 



\ {z 2 + (1 " zf) . (2) 



9 3 z ' q9 ^' 2 

Accordingly, the joint probability density of obtaining a quark or anti-quark coming 
from subsequent decays v — > v + g and g — > q + q at some fixed low Q\ is 

The value of Q v , as dictated by the valon model of the nucleon, is about Q v — 1 GeV. 



For definiteness we take Q v = 0.8 GeV as in Ref. fllll , |i2|1 , which is large enough to 



allow for a perturbative evaluation of the qq pair production. N is a normalization 
constant whose value depends on the flavor of the quark and anti-quarks produced in 
the gqq vertex. 

Once q and q are created, they may subsequently interact with the valons so as to 
form a most energetically favored meson-baryon bound state. The rearrangement of 
such five-component nucleon configuration into a meson-baryon bound state must be 
evaluated by means of effective methods. This is necessary because the interactions 
involved in such a process are within the confinement region of QCD. Therefore, 
non-perturbative interactions take place. Assuming that the in-proton meson and 
baryon formation arise from mechanisms similar to those at work in the production 
of real hadrons, we evaluate the in-proton pion probability density using a well-known 
recombination model approach |14| . 



Within this scheme, the pion probability density in the \ttB) fluctuation of the 
proton is given by 

pm*) = r — r —f(v, aw*, v. *), w 

Jo y Jo z 

where R(x, y, z) is the recombination function associated with the pion formation, 

J j toll) = H,( 1 _v±£), (5) 

and F(y, z) is the valon-quark distribution function given by 

F(y, z) = f3yv(y) zq(z)(l - y - z) a . (6) 

The exponent a in eq. (|5]) is fixed by the requirement that the pion and the 
baryon in the \ttB) fluctuation have the same velocity, thus favoring the formation of 
the meson-baryon bound state. With the above constraint we obtain a = 12.9 and 
a = 18 for the \n + n) and the |vrA) fluctuations of the proton respectively. 

Note that in the original version of the recombination model this exponent was 



fixed to 1 [|14j| . This is basically because in a collision, the only relevant kinematical 



4 



correlation in the model between the initial and final states is momentum conser- 
vation. On the other hand, in the present case the recombining quarks are more 
correlated as they are making part of a single object from the outset. Firstly, meson 
and baryon must exhaust the momentum of the proton [], and secondly, they must be 
correlated in velocity as a bound-state is expected to be formed. 

The overall normalization Nf3a of the probability density P^b must be fixed by 
comparison with experimental data. 

The non-perturbative u and d distributions can now be computed by means of 
the two-level convolution formulas 



d NP (x,Ql) = f 

J X 



dy 



x y 



o 



U-,Ql) (7) 
y 



u NP 



x,Ql) = f ld ^ l -P^{ y )u n {-,Ql), (8) 
Jx y i y 

where the sources d n (x, Ql) and u^ipc, Ql) are the valence quark probability densities 
in the pion at the low Ql scale. In eq. (0), we have summed the contributions of the 
\-K + n) and |7r + A°) fluctuations to obtain the total non-perturbative d distribution. For 
the non-perturbative u distribution of eq. (H), the only contribution originates from 
the |7r~A ++ ) fluctuation. Contributions arising from fluctuations containing neutral 
mesons as the 7r°, are strongly supressed in this model and will not be considered?]. 
We also neglect higher order Fock components. 

The factors | and \ in front of P^a in eqs. (0) and fl8|) are the (squared) Clebsh- 
1 We fulfill this requirement by assuming P^b{ x ) = Pbw(1 ~ x )- See Ref. |li| for a discussion 
about this point. 

2 This supression is associated with their flavor structure (in terms of their parton components). 
In particular for 7r° ~ (dd — uu) fluctuations. In our model the first component arises from the 
recombination of a d quark and a d valon, and the second comes from the recombination of a u 
quark and a u valon. Due to their unflavored structure, neutral in-nucleon quark-antiquark objects 
would tend to be annihilated much more rapidly than charged, flavored ones like ud or du. In other 
words, unflavored qq recombination would be highly inhibited and neutral pion fluctuations would be 
ephemeral. Thus, ir° configurations would be unlikely to occur in comparison with fluctuations, 
although they belong to the same mass multiplet. 
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Gordan (CG) coefficients needed to account for the ~ isospin constraint on the fluc- 
tuation. The CG coefficient corresponding to the \Tt + n) fluctuation is hidden in the 
global normalization of the state. 

We will now compare our results with the experimental data. As the E866/NuSea 
Collaboration measures the ratio d/u at Q = 7.35 GeV, we first compute this quantity 
by means of 

d(x,Q 2 ) = d NP (x : Q 2 ) + q p (x,Q 2 ) 

u(x,Q 2 ) u NP (x,Q 2 ) + q p (x,Q 2 ) • 1 ' 

Here d NP (x, Q 2 ) and u NP (x, Q 2 ) are given by eqs.(^) and (§) and q p (x, Q 2 ) represents 
the perturbative part of the up and down sea of the proton, which we assume to be 



equal. This assumption is exact up to at least 1% |15 



Regarding the difference d — u, instead of computing it directly by subtracting 
eqs. (0) and (H), we will extract it from the d/u ratio as done in Ref. ||. In its 
paper, the E866/NuSea Collaboration employed the following identity to obtain the 
difference: 

d(x) - u(x) = ~ \ [u{x) + d(x)\ . (10) 

d[x)/u[x) + 1 

While the ratio d(x)/u(x) is a direct measurement of E866 the sum u(x) + d(x) 
appearing in eq. (|T0|) is taken from the CTEQ4M parametrization |16[ . 

In Fig. [I], our predictions of d/u and d — u are compared with the experimental 
data from Ref. 0. The curves were obtained using the pion valence distributions of 
Ref. flTj in eqs. (^) and (§) and the proton sea quark distributions of Ref. |1| in 
eq. ©. 

Note that a rigorous comparison of our prediction with the experimental data 
would require that the non-perturbative u and d distributions be evolved up to 
Q = 7.35 GeV. Instead of performing a full QCD evolution program, we pseudo- 
evolve the u NP and d NP distributions by multiplying them by the ratio q(x, Q 2 = 7.35 2 
GeV 2 )/g(x, Q 2 ). The function q represents the corresponding valence quark distribu- 
tion in the proton at the E866/NuSea and the valon scales respectively. This simple 
procedure is satisfactory enough to give us a feeling of the effect of the evolution of 
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the non-perturbative distributions on d/u and d — u 

As can be seen in Fig. [I], the results of the model are impressive, considering that 
we are representing both the difference and the ratio. Nevertheless, in the small- x 
region the model seems to overestimate the value of d — u due to the steep growth 
of the valence quark distribution of the pion as x —>■ 0. If for instance, we multiply 
the quark distribution of the pion by a power of x,f\ the signaled excessive growth at 
very small x is corrected while the rest of the curve does not appreciably changes. 
The d — u difference predicted by the model at the valon scale Q v thus presents an 
inflection point about x ~ 0.05 and goes to zero with x. The description of the d — u 
data is thus improved at the price of having modified the low x behavior of the valence 
quark distributions in pions, a region where they are not well known. In addition, we 
also get a more accurate description of the d/u data in all the measured region (see 
Fig. HD. It should be noted that similar results are obtained by using the low Q 2 pion 
valence quark distributions of Ref. [pOR , calculated with a Monte Carlo based model. 

It is instructive to look at the integrals of the non-perturbative u and d distribu- 
tions in order to get an idea of the relative weights of the |7riV) and |7rA) fluctuations 
in the model. By fixing the normalization of the bound states to fit the experimental 
data, for the unevolved curves in Fig. || (Fig. H) we have Jq 1 dx u NP (x) ~ 0.28 (0.15) 
and fidx d NP (x) ~ 0.47 (0.29). Accordingly, the value of ^ dx [d NP (x) - u NP (x)] 
predicted by the model is 0.19 (0.14) []. This is in good agreement with the experi- 
mental result 0.147 ± 0.039, as given by the NMC fTJ. 

If, on the other hand, we consider the definition of d(x) — u(x) as given by eq. (|T0D, 
our prediction of Jq dx [d NP (x) — u NP (x)] is 0.091 (0.083), in close agreement with 
0.10 ± 0.018, obtained by the E866/NuSea Collaboration |J. Note that this value 
of the integral is significantly lower than the previous one, which we obtained by 
3 A similar strategy has been adopted in Ref.[]l9f 

For simplicity we used a power 1, but other values close to this can do the job as well. 
5 Notice that, as an integral of a non-singlet quantity, L dx [d(x) — u(x)] is independent of Q 2 ||. 

Then, our results at the valon scale remain unchanged after QCD evolution. 
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direct integration of the difference between eqs. (|7|) and (||). This discrepancy is due 
to the modulation introduced by the CTEQ4M u{x) + d(x) distribution used by the 
E866/NuSea Collaboration to extract the d — u distribution Q. 

A similar analysis of the E866/NuSea data has been recently performed in the 
framework of a light cone form factor version of the pion cloud model [^TJ. Predictions 
of this version of the PCM are however not very close to the data. One reason 
may be the use of unnatural hard pion distributions in |7riV) and |7rA) fluctuations, 
which produce large contributions to the u and d distributions beyond x ~ 0.25. 
This drawback in the prediction of d — u translates into the growing behavior of 
the resulting d/u ratio. To obtain an improved description of both d — u and d/u 
within this approach, the addition of an ad-hoc parametrization of the Pauli exclusion 



principle is needed. In particular, in Ref . , the Pauli effect is normalized to 7% while 



the total pion cloud contribution to just 5%. It means that the Pauli contribution 
would amount to a 58% of the total asymmetry. This is a major contrast between 
this approach and the present work. 

Summarizing, we have shown that, including perturbative and effective degrees 
of freedom in a recombination scheme, a pion cloud model alone closely describes 
the recent data of the E866/NuSea Collaboration. With just two parameters, the 
normalization of the |7riV) and |7rA) fluctuations, we have presented an accurate 
prediction of the flavor asymmetry in the light nucleon sea. Remarkably, our model 
results allow an excellent fit of both distributions, difference and ratio, in a consistent 
way. Finally, we have also signaled a possible reason for the apparent discrepancy 
between E866/NuSea and NMC results on the GSR violation. 

Note added: After the conclusion of this paper another PCM evaluation of both 
difference and ratio has been performed Jffi] . In contrast to ours it is based on the 
use of form factors. For the chosen parameters a reasonable fit of the difference is 



6 See also Ref. || for an additional discussion about the discrepancies between E866/NuSea and 
NMC results. 
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obtained but the model predictions for the ratio d/u do not fit the experimental data. 
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Figure 1: Predictions of the model compared with experimental data from Ref. ||. 
d/u ratio (upper) and d—u asymmetry (lower) at Q = 7.35 GeV. Curves are calculated 
with unevolved u NP and d NP distributions (full line) and with pseudo-evolved non- 
perturbative distributions (dashed line) 
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